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In presence of a magnetic field, photons can mix with any particle having a two-photon vertex. 
In theories with large compact extra-dimensions, there exists a hierachy of massive Kaluza-Klein 
gravitons that couple to any photon entering a magnetic field. We study this mixing and show that, 
in comparison with the four dimensional situation where the photon couples only to the massless 
. . . graviton, the oscillation effect may be enhanced due to the existence of a large number of Kaluza- 

, Klein modes. We give the conditions for such an enhancement and then investigate the cosmological 

■ and astrophysical consequences of this phenomenon; we also discuss some laboratory experiments. 

' Axions also couple to photons in the same way; we discuss the effect of the existence of bulk axions 

in universes with large extra-dimensions. The results can also be applied to neutrino physics with 



(N 
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^ ' extra-dimensions 



I. INTRODUCTION 

> . 

' It is well known |jl|,^ that photons can be converted into gravitons by a magnetic field in a standard four dimensional 
^ ^ spacetime. The propagation eigenstates are then mixtures of photon and graviton interaction eigenstates. In quantum 
words, this mixing is due to the fact that any particle which has a two-photon vertex can be created by a photon 
[ entering an external electromagnetic field |^,^ and can oscillate coherently with the photon. Classically, this can 
I be understood by the fact that an electromagnetic plane wave cannot radiate gravitationally in vacuum since its 
' stress-energy tensor contains no quadrupole H. But, a time varying quadrupole appears (due to interference) when 
an electromagnetic plane wave propagates through a constant magnetic field The implications of this effect 

O ^ on the cosmic microwave background (photons) has been considered and it has been shown that it will be 
I undetectable for standard cosmological magnetic fields A similar effect also generically happens for axions (and 
for any particle having a two-photon vertex). The pho ton-axion (see ]9|-pT|| for reviews on axions) mixing has yet been 



X 



studied in details by many authors (see e.g. 12|- 14]) and is used in experiments, since the pioneer work by Sikivie 



to put constraints on the axion parameters |15|-18[| (see e.g. for an up to date review on such experiments). 
Recently, a lot of interest has been raised by models where the universe has large extra-dimensions ||2^,|2l[ . In such 
models, the Planck scale, M4, is no longer a fundamental scale but is related to the fundamental mass scale of the D 



■ dimensional theory, Mjj, through |21 



M| = R''M^+^, (1) 



where i? is a length scale (usually taken to be the radius of the n = D — 4 compact extra-dimensions) and M4 = 
M4/-\/87r. Md can be significantly smaller than M4 at the price of having large extra-dimensions. These ideas can be 
naturally embedded in fundamental string theories with a low string scale p^-|2^ (see also [|6| for earlier discussions 
on TeV scale extra-dimensions). In these models, gravity can propagate in the D dimensional spacetime (bulk space 
time) whereas the standard model fields are localised on a 3-brane. An effect of the compact extra-dimensions arises 
from interactions between the Kaluza-Klein (KK) excitations of the gravitons (or other bulk fields) which are seen 



in four dimensions as a tower of massive particles 28 1. Constraints on the size of these extra-dimensions can be 
obtained both from the laboratory physics and from astrophysics and cosmology (see e.g. |29)). For instance, the 
emission of KK gravitons induces an energy loss in many astrophysical objects such as the Sun, red giants and 
supernovae SN1987Aj|ltl implying the lower bound > 30 - 130 Tev (2.1 - 9.2) TeV for the case of n 2 (3) 
extra-dimensions [^J31| . Some authors |20 2^p^,Q also have recently investigated the possible presence of axions 
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in the bulk which would be coupled to the brane degrees of freedom. Such a bulk axion also gives rise to a tower of 
KK-states as seen from a four dimensional point of view. 

The goal of this article is to investigate the effects of the photon-KK graviton and photon-KK axion oscillations 
and to estimate their effects in cosmology and astrophysics, as well as terrestrial experiments. Since there is a large 
number of KK states with which the photon can mix, one can expect a departure from the usual four dimensional 
result. We first describe (§ |ll|) the photon-KK graviton system starting from a D dimensional action, linearising it 
and compactifying it to four dimensions. We also show (§ III) that the photon-axion mixing is described by the 
same formalism and lead to the same effects as for gravitons. Then, we turn to investigate the mixing in itself. For 
th at pu rpose, we sum up the known results of the mixing of a photon with a low m ass particle in four dimensions 
(§ IV A ) and then discuss the most general case of a D dimensional spacetime (§ [V B| ). The general expression for the 
oscillation probability is then evaluated in the particular cases of a five (§ ^) and of a six (§ [V^ ) dimensional spacetime. 
We show that, as long as one is in a weak coupling regime, one can add the individual probabilities which leads to an 
enhancement of the oscillation probability if the KK modes are light enough. We also show that there exists a regime 
where the photo n mix es strongly preferentially with a given KK mode. We generalise our results to an inhomogeneous 
magnetic field (§ VII) when the probability of oscillation is small, and then turn to the cosmological and as trophys ical 
situations where such effects may be ob served . We study the case of the cosmic microwave background (§ VIII A), of 
pulsars (§ VIII B ) and of magnetars (§ VIII C| ). We show that even if the enhancement of the oscillation probability 
can be very important, it is still very difhcult to observe this effect in known astrophysical systems. To finish (§ IX) 
we discuss laboratory experiments and particularly polarisation experiments. For that purpose, we describe the 
computation of the phase shift between the two polarisations of an electromagnetic wave and compare the result to 
the standard four dimensional case. 



II. EQUATIONS OF MOTION FOR THE PHOTON-GRAVITON SYSTEM 



Following |2^,|2^ , we consider a field theory defined by the D dimensional action 

1 



Sr 



where gAB is the D dimensional metric with signature (— , +, +), k^^ = SttGd = Mj 
Lagrangian. The indices A, B, ... take the value 0, ..3, 5.., D and we decompose as 

z^ = ix^,y'^) 



-{2+n) 



(2) 

and Cm is the matter 
(3) 



with ^,v, . . . = 0, . . . , 3 and a, 6, . . . = 5, . . . , Z?. 

This theory is considered as being a low energy effective theory valid below some cut-off A/max in energy (see e.g. 
p0| , pl|j25| ) . We will discuss in this paper only the cases n = l and n = 2 in details, and our conclusions, in these cases, 
are mostly cut-off independent. We stress that the relationship between this cut-off and the fundamental string scale 
(if one wishes to embed these theories in superstring models) can be much more complicated than what is naively 
expected (see Q). 

We expand the metric around the D dimensional Minkowski spacetime as 



h 



gAB = riAB 



AB 



M 



l+n/2 
D 



(4) 



where rjAB is the D dimensional Minkowski metric. Inserting (^ in (^ and using the definition of the stress-energy 
tensor as ^J^Tab = 2S{Cm\/--g)/Sg^^ , so that Cm\f--g = Cmo — h^^TAB /2Mj^^^^^ , we obtain the linearised action 



S. 



D 



d'^Z 



ih^^d'^dchAB - \hid^dchl + ^h^^'dAdBh^c + \hidcdBh'''' 



h^'^'dAdch^ 



1 



l+n/2 



M 



D 



(5) 



We compactify this theory to get a four dimensional theory and use the periodicity on y°' to expand the field h^^ as 



hAB{z^)^Y.'^^'' 



exp [i- 



Vn 



R 



(6) 
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where Vn = (27ri?)" is the volume of the compact n dimensional space (assumed to be a cubic n-torus). Hab is split 
into a sum of KK modes living in the four dimensional spacetime. The ordinary matter being confined to the brane, 
its stress-energy tensor must satisfy 

rAs(;^^) = T^.(x^)<S("Hy^)«- (7) 

In which follows, we restrict our attention to the case of an electromagnetic field Ff^i, for which the stress-energy 
tensor is given by 

T^. = F^xF,^ - jV^^F^PFxp. (8) 

The fields 

^^AB '^^^ be decomposed into spin-2, spin-1 and spin-0 four dimensional fields |^,^. Only spin-2 
and spin-0 particles couple to ordinary matter and spin-0 particles couple only to T^. For the electromagnetic field 
= classically so that the only relevant KK modes (at the tree level analysis of this article) will be the spin-2 
particles for which the action (||) reduces to ]2^ , ^ 



5*4 = d 



where m?. — / is the square mass of the KK graviton, □ = d^d^. 

The equations of motion deduced from the Lagrangian (||) are the coupled Einstein-Maxwell equations 



(9) 



(□ - ml)G^$ = (10) 

d^Gf. = (11) 
G(f ^ = (12) 
2 

M4 



d^F"'^ {g^p^°"'F^ - G^P^^-'F^") = 0. (13) 



When p 7^ 0, the field g'^^ has 10-1-4=5 de grees of freedom which is what is expected for a massive spin~2 particle. 

We now consider an electromagnetic plane wave in the presence of a magnetic field Hq which is assumed constant 
on a characteristic scale Ac in the sense that its variation in space and time are negligible on scales comparable to the 
photon wavelength and period. We define the basis 

k Ho± 

such that {cy , e+ , e|| ) is a direct orthonormal basis of the three dimensional space and where Hq± is the perpendicular 
component of Hq with respect to the direction of propagation k. We decompose the KK gravitons in scalar (S), vector 
(V) and tensor (T) as 



(V) G^^v}^^vf4 + V^e^, Gg^^O, Gf = ^k^,V^^ (16) 
(T) Gi^^O, G®=0, Gf=Gfe++Gfe^, (17) 

where = (jtikj — ^Sij^, i,j = 1..3 and a dot refers to a time derivative. The polarisation tensor of the graviton 
modes efj is defined by 

e^_,. ^ (efe- - e+e+) 6^, + 2ei+efsi. (18) 
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The advantage of such a decomposition is that the scalar, vector and tensor contributions decouple. The five degrees 
of freedom of each massive spin-2 KK gravitons have been decomposed in one scalar mode (0^^-*), two vector modes 
(K^/^x) ^^"^ tensor modes (G^j^). Each of these modes satisfies independently the constraints (|lTl-|l2|). 
We consider an electromagnetic wave with a potential vector of the form 

A^i{A^{u),A+{u),{))e-'^\ (19) 

where u is the coordinate along the direction of propagation. We have introduced the arbitrary phase i so that its 
electric and magnetic fields are 



E = -dfl = {LuAy, {u),ujA+{u), 0)e-*'^* 

B = cm\{A) = {^iduA+{u),iduAy,{u),0)e-''^K 



(20) 
(21) 



The stress-energy tensor of these waves in the presence of Hq has no vector component. Its tensor component is given 
by 



A=+,x 



(22) 



We see, as expected, that a plane wave possesses a tensor part only if it propagates in an external field and that 
the polarisations + and x of the electromagnetic wave couple respectively to the polarisations + and x of the 
gravitons. The electromagnetic wave generates also scalar perturbations, but it can be shown [|j that (in the usual 
four dimensional case) the total energy converted in this scalar wave are negligible compared to the tensor contribution. 
In the following, we concentrate on the tensor modes. 

The equation of evolution of this system is given by the Einstein equation (|lO|) which reduces to 



and the Maxwell equation (O) which reduces to 



M4 



■duAx, 



(23) 



(24) 



where we have used the ansatz Gij = J2\ G\{u)e' 



ij 



for the gravitons. 



Since we have assumed that the magnetic field varies in space on scales much larger than the photon wavelength, 
we can perform the expansion co'^ + = (ui + idu){u} — idy) — {lo + k){Lo — i9„) for a field propagating in the +u 
direction. If we assume a general dispersion equation of the form lo — nk and that the refractive index n satisfies 
\n — 1| <C 1, we may approximate lj + k — 2uj and k/uj — 1. This approximation can be understood as a WKB limit 
where we set A{u) = I4('u)|e''''" and assume that the amplitude \A\ varies slowly, i.e. that du\A\ <C k\A\. In that 
limit, the system (P3|-p4|) reduces to 



the matrix M\ being given by 



r Ax 



- idu+Mx] 



G 



(0) 



G 



(9) 



Aa 


Aa/ 


Aa/ 




Aa/ 


Aa/ 












Aa/ 





^(1) 








Aa/ 










^(9) 



(25) 



(26) 
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with 



and A(^)^p^^)A„. (27) 



M4 



is a n-uplets (p^^ , Pq^ , . .pq^ ) of integers and we have ordered the Am such that 

a(9) < a(9+1) 

and A„i is defined by 

Each Am' appears witli a multiplicity given by the number of n-uplets having the same norm J2i=i nPqi- We define 
the two series {ri)i>i and (si),;>i such that 

ACn-i) ^ A (n) _ A(n+i) _ A (n+s,-!) A (n+s,) _ A (n+i) /r,qN 

(r) 

We have n+i = j'i + s^, and Si is the multiplicity of the element Am , i-e. the number of times it appears in the 
matrix (|2^). ri is the rank in the series (Am\ Am\ ■ • ■) where the i**^ distinct value of Am^ appears for the first 
time. In the case of a five dimensional spacetime one can easily find out that si = 1 and Si = 2 for i > 1 and that 
ri = 0, r2 = 1, ra = 3, . . . In the case of a six dimensional spacetime, Si = (1, 4, 4, . . .), — (0, 1, 5,9,.. .). Introducing 
the cut-off Minax discussed above, we require m'i = jP' /R^ < M^^^, which using (|]) translates into < _p^ax "^i^h 

Setting Mmax Md, we obtain Pmax ~ (^Mi/MD)^^" . For n = 2 and Md ^ ITeV one finds 

Pmax - 10^^ (31) 

which means that we have to consider a very large number of KK states. We also define a maximum index, N say, 
for the series Aif defined by 

N = sup{q I P(q)=pLx}> (32) 

which translates into a maximum index Nj:, for the series Si and . We stress that the number of KK modes relevant 
for the photon- KK graviton oscillation is likely to be smaller than TV due to decoherence effects, such as the source 
and detector finite width in momentum, the wave packet separation for massive (and non-relativistic) KK modes... 
(see e.g. |Q for a description of these effects in the case of neutrino oscillation). 

The term A a can be decomposed as A a = Aqed + Aqm + Apiasma- The first term contains the effect of vacuum 
polarisation giving a refractive index to the photon (see e.g. Adler |^) and can be computed by adding the Euler- 
Heisenberg effective Lagrangian which is the lowest order term of the non-linearity of the Maxwell equations in vacuum 
(see e.g. [^^) to the action (^) ^. The second term describes the Cotton-Mouton effect, i.e. the birefringence of 
gases and liquids in presence of a magnetic field and the third term the effect of the plasma (since, in general, the 
photon does not propagate in vacuum) . Their explicit expressions are given by 



^plasma 



plasma 



2uj 



A^^ - A+ = 2TTCHI (33) 



^The equation of motion derived from ^ is ( ]2^ ) with Aa = 0. We intentionaly omit the Euler-Heisenberg contribution in 
the presentation for the sake of clarity. Its Lagrangian is explicitely given by Ceh = \_{F'^'^ F^u)^ + ^(^"'"'-Fmi^)^] . 
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with ^ = {a / A5'k){H[)± / H c)"^ , He = rn1/e = 4.41 x 10^"^ G, me the electron mass, e the electron charge and a the fine 
structure constant. C is the Cotton-Mouton constant |Q; this effect gives only the difference of the refractive indices 
a nd th e exact value of C is hard to determine |39] ; we will neglect this effect but for the polarisation experiments (see 



[X B ) . The plasma frequency Wpiasma is defined by 



^^piasma = 47ra — , (34) 



rie being the electron density. Note that the Am' are always negative whereas is positive if the contribution of 
the vacuum dominates and negative when the plasma term dominates. 

The equation of motion (|2^) reduces to the one studied by Raffelt and Stodolsky |^ when one considers only four 
dimensions so that A4\ contains only the massless graviton and is then a 2 x 2 matrix. The main difference lies in 
the fact that now the electromagnetic component couples to a large numbers of KK gravitons. This can be compared 
to some models of neutrino oscillations in spacetime with extra-dimensions . We should also note that the two 

polarisations are, as expected, completely decoupled and obey the same equation of evolution. In the following of this 
article, we omit the subscript A of the polarisation. 



III. EQUATIONS OF MOTION FOR THE PHOTON-AXION SYSTEM 



Before turning to the study the photon-graviton mixing, we consider the case of axions and show that the photon- 
axion mixing can be described by the same formalism. 

We consider the generic action [p2|j3^ for the bulk axion photon system 



5*4 



E 



fPQ 



1 



(35) 



where the are the mass eigenstates of the axions and mp their masses. F^^y = \e^ypaFP'^ is the dual of the 
electromagnetic tensor, e^^pa being the completely antisymetric tensor such that £0123 = +1- As for the gravitons, 
the mass spectrum is expected to be discrete, the states can be labelled by a n-uplet p and is expected to have a 
typical spacing of 1/R. We have considered here that every axion KK state a'^^ couples to the photon with the same 
coupling l//pQ. This is only expected to be true if the typical mass, mpQ, given to the axion zero mode by instanton 
effects is much lower than the typical KK mass 1/R . Let us further stress here that for such bulk axions the usual 
relationship between the axion mass and the PQ scale does not hold anymore, so that one expects to see interesting 
new effects to appear Q. Inspired by the usual bounds on /pq, we take /pq of order 10^° GeV. However we stress 
that the usual bounds on /pq may be modified partly because of a large number of axion-like particles coupling to the 
standard model fields. For example, we expect that the astrophysical bounds will be more stringent mainly because 
a star will now be able to emit all the energetically accessible modes (see |Q and also |Q for a discussion on relic 
axions oscillations). 

We do not consider the perturbations of the metric and work in Minkowski spacetime since we are interested in 
the interaction between the photon and the axion. We deduce from ( |35|) the coupled Klein-Gordon and Maxwell 
equations 



/PQ 



/PQ 



(36) 



(37) 



We now decompose the electromagnetic wave as in (19-pl|) with respect to the basis (p^), so that the former system 
reads 



[D - mt) a 



AHoi 
/pq 



■A, 



(38) 
(39) 
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where we have decomposed the axions as a^P>{u) exp {—iujt). Using the same WKB hmit as in section IV, we obtain 
the hnearised system 



The matrix A4 is now defined by 



{lu ~ idu +M) 



M - 



Ax 

nip) 



= 0. 



/A+ 














Ax 


Am 


Am 





Am 


^(0) 








Am 








A 19, 



(40) 



(41) 



with 



Am = 



4go± 

/PQ 



A(p)-_!i- 
2u;' 



(42) 



A+ and Ax being given by equation (p3|). This system reduces to the RafFelt and Stodolsky Q system when we 
consider only four dimensions. Only the component x, i.e. parallel to the magnetic field, couples to the axions. This 
is a major difference compared with gravitons for which both polarisations of the photon evolve alike whereas here 
only Ax is affected by the mixing. 

One of the goal of thi s se ction was to set the theoretical framework for further experimental studies of photon-bulk 



axion oscillations (see 



IX 



for a more detailed discussion) and to show that it is described by a similar formalism as 
photon-KK graviton oscillations (under the validity conditions explained below equation (p5|). 



IV. PHOTON-KK STATE MIXING IN A HOMOGENEOUS FIELD 



We now describe the physical implications of the system (|2^-|24|) and start by re viewi ng briefly the well studied 
problem of the mixing of a photon with a low mass particle in four dimensions (§ LVA). We then give the exact 
expression of the oscillation probability in D dimensions and discuss qualitatively its magnitude and the effect of the 
coupling of the photon to a large number of particles. 



A. The usual photon mixing with a low mass particle 



This case was well studied in the literature, see e.g. Raffelt and Stodolsky IJ] and we just summarize the main 
features of the results to compare to the case of a spacetime with extra-dimensions. For the mixing with a single 
particle of mass m, the matrix M reduces to 

with A„i = —rv?j2L0. The solution to the equation of motion ( ^ ) is obtained by diagonalising Al throught a rotation 

"A' 
G 

with 



cos I? sin -Q 
— sin cos 





A 


) 


G 



(44) 
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and where we have introduced a = A^/A^ and /3 = Am/A>,. We obtain by solving ( |25| ) in this new basis 

=e~'^^"^'(0) 

G'{u) = e-'^>G'{Q) (46) 
where a global phase uju has been omitted. The two eigenvalues A^ and A^ of are explicitly given by 

Going back to the initial basis, we obtain 

A{u) = (e-^'^A" cos^ 7? + e-'^> sin^ i?) A(0) + sin i? cos d (e"*^^" - e-^^>) G(0), 

G(m) = sini^cost? (e""^'^" -e-*'^>) ^(0) + (e^'^^s" cos^ + e"*'^'^" sin^ G'(O). (48) 

The oscillation probability of a photon into a graviton is computed by considering the initial state (^(0) = 1, G(0) 0) 
and is given by 

P(7 ^ g) =1 (A(0) I G{u)) |2 = sin2 {2d) sin^ f , (49) 



with 



2 

(Aosc^VS) 



2 sin^(Ao,cM/2) 



A..e ^ A^ - A' = = ^ = I^Aa (51) 

^ 3 cos2i9 sin2z? cos2i? ^ ^ ^ 



The oscillation length is thus given by ^osc = 27r/Aosc- We see that the oscillation probability cannot be larger than 
(Aj\/u)^. One has to be aware that d depends on the polarisation index A. 
It is interesting to single out the two following limiting regimes: 

• The weak mixing regime in which ^ 1 so that the probability (^) reduces to 

^^(7 - 5) = 4^^^ sin^ (^^^^") ■ (52) 

When the oscillation length l^sc = §^ is large with respect to the coherent path distance u, the weak mixing 
probability can be further approximated (with Amu <C i9 <C 1) by 

P(7 ^ g) ~ [aAxuf ~ (Amu)' • (53) 

• The strong mixing regime in which the mixing is maximal, i.e. when ~ 7r/4, so that the oscillation probability 
reduces to 

P(7 -> g) = sin^ (Amu) (54) 

and the oscillation length to 

^osc — "T • (55) 

A complete transition between a photon and the light particle is then possible. This can only happen when Am 
and A A have the same sign (see equation (^5|)). We further note here that the width in /3 of the strong mixing 
region is roughly given by a according to equation (Ha). 
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B. Mixing in D dimensions 



1. General result 



To compute the oscillation probability in a spacetime with extra-dimensions, we first have to solve ( p5| ) which implies 
the diagonalisation of the matrix (p6|). We present the explicit and detailed computation of both the eigenvalues and 
eigenvecto rs i n appendix |^. We then compute in appendix ^ the explicit form of the oscillation probability (see 
equation (B5)) 



P(7 ^g)^l 



Nd 

E 

i=l 



Taking into account the fact that /^^^ = 1, it can be rewritten as 



iVr 



where the coefficients f'^. are defined by (see equation (p6|)) 



(56) 



(57) 



f 



Ne 



~[ iVj - Pi) 



(58) 



The expressions (|^ and ( psf ) depend on the eigenvalues yi solutions of the equation (|A9|). Introducing the notations 
y = x/Aa, a = Am/Aa, (3 = Am/A\ and Pi = A^'V^a, the eigenvalues equation (|A^) can be rewritten as 



No 



(59) 



The photon-KK graviton oscillations is then completely described by the set of equations (|57|-|5S|). 

Indeed, it is difficult (even impossible if n > 1) to compute analytically the roots of (p9|). For instance the coefficients 
Si are not known analytically if n > 1; it is of course possible to compute P(7 g) numerically, but this is not our 
purpose here. In the next two sections we derive the oscillation probability in the two cases n = 1 and n = 2 in a 
range of parameters dictated by the systems where such a mixing may appear. In the next paragraph, we discuss 
qualitatively the results found there, stressing some new effects due to the presence of a large number of KK states, 
as well as to the degeneracy of each KK level. 



2. Qualitative discussion 



We only discuss the cases where a is small in comparison to (3 as dictated by the physical systems studied in § VIII 



and § [X. Two different limiting regimes appear, a large radius regime (when \f3\ is smaller than unity) and for which 
there is a significant effect of the extra-dimensions, and a small radius regime (when |/3| is larger than unity) and for 
which there is only small departure from the usual four dimensional photon mixing. 

Let us first discuss the large radius regime where is smaller than unity. As in four dimensions, according to the 
respective value of Ax and of the Am'' 's, two behaviours can appear: 

• A strong mixing regime with one given KK state , if there exists a state K such that 

|AW_Aa|« V^Aa/. (60) 

This can happen only if /3 > 0, i.e. when plasma effects dominate over the vacuum polarisation. We stress also 
that since we have assumed along this discussion that a is lower than (3 there is at most one KK state which 
can mix strongly with the photon. The total probability will be found to be dominated by a term of the form 
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^^(7 ^ 5) = (1 - ^) sin^ 



Sja^ . , / 1 - ft 



4 2. 1 ^^^^ ) ' (61) 



with 

Aifj = yi]?AM (62) 

(see § ^ and § ^ for a detailed derivation), rj is much smaller than unity. As will be shown later, this form 
accounts for keeping only the dominant part of each /J. . Depending on the argument of the sines, the probability 

is either dominated by sin^ Ao^''u or by the correction coming from the modes i ^ K. This shows a first 

departure to the four dimensional case due to the degeneracy of the KK level K\ the oscillation length associated 
with the strong mixing state (labelled by K) is lowered by a factor ^sk which can be very large. Moreover, the 
width of the region in Aa of strong mixing is larger by a factor ^sk than in the usual case [see below equation 

(H)]. 

An other important difference with the usual four dimensional situation, where the strong mixing regime can 
only occur when A^ crosses the unique A^ characteristic of the mixing state, we now have more possibilities to 
be in that regime, where a complete transition between the photon and the graviton is possible. Because of the 
presence of a KK state Am^ in any interval in A^ of typical width A,„, only fluctuations of A a of order A™ can 
lead to it. 

• A weak mixing regime where for all g, — 1| 3> a. The oscillation probability is then given by 

This contribution is exactly the one that will be intuitively thought of and obtained by summing the individual 
oscillation probabilities (p3) of the photon into each KK state with the mixing angle 



2a 

tan2i?, EE — — -. (64) 

-l^ Pq 



There are roughly three contributions to the sum (|63| ) that we estimate as follows. 

1. All the states such that \[3q\ <^ 1 mix with the photon with approximatively the same angle dq aii we 
neglect Pq with respect to unity in (|4|). The order of magnitude of the probability of oscillations with 
these states is then 

,2„- 2 ( Aa^ 



P(7 -> g) ^ ANia^ sin^ yY"") ' ^^^^ 

Ml can be estimated by counting the number of modes such that Pq < P^i — 1 with their multiplicity (see 
appendix i.e. Mi - Xlf^i ^""^ ~ ^i" ~ /^""''^ so that 

P(7-3)-^Bin^(^«). (66) 

This already shows that the oscillation probability can be greatly enhanced (by a factor /3~"/^ with respect 
to the four dimensional case with the same mixing parameters [obtained from equation ( p2[ ) with |/3| <C 1]). 

2. All Pq such that \Pq\ ^ 1 have a mixing angle roughtly estimated by i?g ~ ct/ Pq and their contribution to 
the probability is of order 

P(7 ^g)^Aa'Y.^2 ■ (67) 

This series is difficult to evaluate since the oscillation length is different for each KK state. When n < 3 it 
can be bounded by a^//3"/^ so that this contribution is at most of the same order of magnitude than the 
previous one. 
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3. The contribution of the (3q such that Pq ^ 1 which only exists if /3 > is bounded by Sft^i -^i/iPi ~ 1)^- 
First of all, since |/3g — 1| ;:g> we are never in a strong mixing regime. Now, we single out (3k, the closest 
fii to unity from which it follows that Vi ^ if, \(3i — 1| > 13/2 and thus the contribution of all the /3i ~ 1 for 
i 7^ if is bounded by {a^ / 0^) X]/3 ~i i^K '^^'^ dominated by the contribution of the term K given 

by a'^SK/\l3K — Ip according to the relative value of \(3k — 1| in imits of (3. 

In conclusion, the weak mixing case is characterised by an enhancement of the probability by a factor at least 
p-n/2 |-j^g ^j^g £g^g^ that the photon couples to a large number of KK states. We further note that when /3 < 
and <C 1 one can obtain an absolute bound on the oscillation probability of order a^//3^ for n < 3 (when 
n = 2, this bound is given by lOQaV/32 [see appendix @]). 

We now turn to the small radius regime where 3> 1 and in which the photon mixes preferentially with the zero 
mode. The probability (^) can be expressed as 

^^(7 - 9) = (1 - e)PMl - 9) + 4E (^^^ sin^ I^ ^a(I-A) ^^ (gg) 

where P^d is the oscillation probability for the mixing with the zero mode and is given by ( |49| ) and e ^ 1 is the 
correction of the oscillation probability with this mode coming from the existence of the extra-dimensions. In this 
case, the lightest massive KK mode is so heavy compared to the photon effective mass that it can barely be excited by 
the photon. The contribution of the other KK modes can be shown to be bounded by 0(0^ The contribution 
of the massive KK states is suppressed by a factor 3> 1. 

Introducing the Compton wavelength, say, associated with the effective mass of the photon and defined by 

= \cjA^\-^/\ (69) 

the required condition to be in a large radius regime can be rephrased as A^ < R, i.e. that the average scale associated 
with the photon is smaller than the radius of the extra-dimensions. The latter is expected to be of the order of the 
centimeter for two extra-dimensions. 

In the two following sections, we derive these results in details for a five and six dimensional spacetime. Let us 
stress here than when n > 3 the oscillation probability will strongly depend on the cut-off in energy in which case a 
more precise knowledge of the whole theory and the exact experimental situation (to take into account decoherence 
effect) are needed. We emphasize that in the following we have set the cut-off to its maximum value in order to be 
very general. The computed probability is thus the maximum one and the bounds on the parameters used to derive 
it are the most stringent. 



V. ESTIMATION OF THE PROBABILITY IN A FIVE DIMENSIONAL SPACETIME 



In this section, we present the computation of the eigenvalues and of the oscillation probability when n = 1. As it 
will be seen the computation is easier in this case because the sums in ([5^p9|) can be expressed in terms of circular or 
hyperbolic functions. Although the case n = 1 is generally regarded as in contradiction with observation (see however 
e.g. p^), this computation will teach us a lot about the mixing with a large number of particles. 

In a five dimensional world, the mixing matrix A4 is cxplicitely given by 



Mx = 



/ A, 
Am 
Am 
Am 

AAf 



Am 









Am 



-1 








Am 



-1 




Am 




-4 





Am 





-4 

2WZ 



\ 



\ 



(70) 



This matrix can be compared to the one obtained for neutrino oscillations is spacetime with extra-dimensions (see 

is twice degenerated so that ri = 1 , = 2i — 1 and Si = 2 



e.g. IQ). We see on that example that for 5 7^ each A 
for i> 1. 



(«) 
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The characteristic eigenvalues equation det(A^ — xl) = 0, in that simple case, can be obtained by developing the 
determinant of order 2N + 2 with respect to the last line to get a recursion relation with the dete rmin ant of order 
2 A'' and to find the limit of this series. Indeed it leads to the same result that the general equation (A5). 

With the notations of the former paragraph, the eigenvalues equation (A5) can now be rewritten after resummation 
(see 1.217 in ||) as 



2/-l = -/C(y//3), 



where the function K. is defined by 



'^y'^'-'^m \ _coth^|a;|i/2 (2;<0). 



(71) 



(72) 



The oscillation probability is then given by (|57|) where the coefficients (pSh are now reexpressed after resummation as 



1 -1 



where the function X = X^fcezl^ ^ ^ obtained from ^73) as 



I{x) 



cot ^T^/x + 7ry/x{l + COt^ TTy/x) 

2|x|3/2 \ coth^T^/\^ ~ n^/\x\{l - coth^ n^/\x\) 



(x > 0) 
{x < 0). 



(73) 



(74) 



We now determine an approximation of the solutions of the eigenvalues equation (|7l|) and then of the oscillation 
probability in the cases /? < and /3 > assuming that <C 1. 



Setting /3 = —f3 > 0, we have /3j 
such that 



A. /3 < 

1. Eigenvalues 

-{j — 1)^/3 and one can easily sort out that the solutions of equation (|7l|) are 



2/1 > 1, yj+i 



with 



vi solution of V — 1 = TT^^i — cothvr. 

^ >-■ f 1 + f-y 

?7i>l solutions 01 1 — i/^TT^'l/ COtTTi/^^. 



(75) 

(76) 
(77) 



In figure |^, we depict the graphical resolution of this equation. 

The resolution of equation (|7^) then splits into the three following cases: 

1. Ui'. When yi/P <^ 1, (|76| ) reduces to y~l — /y so that yi ~ 1 + and one can then check that the condition 
yi/P <^ 1 is equivalent to /3 » 1. 

When » 1, we set 2/1 = 1 + e with e > and (^ implies that < e ^ Tra^{l + e)-^!"^ j 0^1"^ < na'^/P^^^. 
Then if ^ 1 

we deduce that yi ~ 1 + -ko? j 0^1"^ and that the initial condition on y\ is equivalent to 

^<Cl. 

In conclusion, if c? j 0^1°^ < 1, 



2/1 



l + a2 (/3>1) 
l+TT^ (/3«1). 



Note that these two solutions can be rewritten under the more compact form 

2/i^i + ^/c(r'). 



(78) 



(79) 
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2. j/j+i; j > 1: Since yj+i satisfies 



U-i)<d-f^<j 



we set J —yj+i/ P — {j — 1) + Cj+i with < Cj+i < 1 and equation { \n\) rewrites as 



1 



J -I 



na cotTrej+i 



(80) 



Now, if / P ^ 1, tlie l.h.s. of (|80| ) being larger than unity, it implies that cotncj+i ^ 1 which thus behaves 

2a2 



as l/nej+i. We can then solve (KQ) for ej+i to get 



1 - 

This expansion is valid whatever the magnitude of (3 as long as a^//3 ^ 1. 
3- y2'- y2 is the solution of (|77| ) such that < \/ —yijP < 1/2- Setting z = —y2/(3, (|77| ) leads to 

cot TT-y/z 



(81) 



1 + /3z = TT- 



/3 



(82) 



with < < 1/2. The l.h.s. of (|8^) being greater than unity, it implies that, when a^//3 ^ 1, cot(7r\/z)/Y^ ^ 
1 and thus behaves as I/tt^/z. At lowest order ( ^2[ ) then leads to z ~ a^//3 and then 

,2 



2/2 ^ -a 

Again, this solution is valid whatever (3 such that o? jP <SC_ 1. 
4. Summary: When ^ 1, the roots of ( |7l| ) are well approximated by 

yic.i + ^/c(r') 



(83) 



(84) 



for all /3 such that ^ 





J L 



0.01 

0.1 

0.5 



0.5 1 1.5 



-0.5 

y 

FIG. 1. Example of a graphic resolution of equation (|7l|) where we have chosen o? = (0.01, 0.1, 0.5) and /3 



-0.5. 
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2. Oscillation probability 

Assuming that a/|/3| <C 1, we can now expand (fz^ ) to get the foUowing behaviours of the coefficients 

fl^l-^Iir^), (85) 

(87) 

One can easily check that, at this order, ^ /J. = 1. Using the form (^7|) of the osciUation probabihty we deduce that 

F(7 - 5) ^ E sin^ (^^^^-) ' (88) 

as announced in (p3|). It can be checked that the dominant contribution to the probabihty comes from the terms 
flP., in @. 

It is worth noting that the cut-off of the theory does not enter the resuh, due to the fact that in this special case 
all the sums arc converging. Since Pj < and a ^ 1, we are always in the weak mixing limit and the oscillation 
probability is well approximated by the sum of all the individual oscillation probabilities. The individual oscillation 
lengths are given by 

fU) = ^ < ^ =/i) 



B. /3>0 
1. Eigenvalues 

Since /3j = [j — 1)^/3, one can easily show that the roots of ( |7l| ) are such that 

yi < 0, /3j < %+i < I3j+i (89) 

with 

2/1 solution of 1 — y — tt—\ — cothTrW — — , (90) 

p ^ -y VP 

solutions of w — 1 = TT— W — cot TT. / — . (91) 
P ]] y \l 13 

In figure || we depict the graphic resolution of this equation. We introduce K the index of the closest f3i to unity. 
Contrary to the previous case, the discussion has to be split in four steps: 

1. yi: When —yi/(3 <C 1, ( po| ) implies that 1 — yi ~ —a'^/yi so that yi ~ —a^ and the initial condition reduces to 
q;^//3 <C 1. Thus when a^//3 <C 1, whatever the magnitude of /3, 

2/1 ~ -a^. (92) 

2. yj,l < j < K: We set \/yj/ (3 = {j — 1) — £j with 1 > ej > 0, so that ( ^ rewrites as 

7ra^ cot 7re.; 

Since the l.h.s. of ( |9^ ) is positive we deduce that ej < 1/2. Now, taking into account the fact that 1 e 
]{f^K + /3if-i)/2, (Pk + Pk+i)/21 we deduce that 1 - y^- > 1 - yK-i > {(3k - (3k-^i)/2 = (3{K - 3/2)/2 and 
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thus (since K > 2) that 1 — yj > P/A. If o? j ff^ ^ 1, then ( |9^ ) imphes that cot Trcj ~ l/Ti'tj ^ 1 from which we 
deduce tj and then 

y^<^<K - - 1)' - (94) 



3- j > ^"^^ + 1: The argument foUows the same hnes as the previous one but we now set \/yj / (3 = [j — 2) + €j 
with 1 > ej > 0. We can now deduce from yj — 1 > yK+2 — 1 that yj — I > (3/2 and then that if o? j 0^ ^ 1, 

2/,>K+i ^ Pij 2f + (95) 
4. yx, yx+i'- If > 1, we set y^ — f^xi^ + e) and y^+i = + e')i '^^ can use the property of equation (|9l|) 



to deduce, as before, that 



K-3/4 , „ , is: -5/4 
0<^<(F^ 0<-6'<^^^^-^ (96) 



and then if /3 is smah, we can conclude that e and e' are small compared to unity. Setting 6 = — 1, e and e' 



are solution of (91) which reduces to 



S + e~— =^ 2e~6± + Sa^. (97) 
e 

Now, if (5 > 2V2a, we deduce that 

Pk - I Pk 

depending on the sign of S and with the constraint yx < yx+i- On the other hand, if (5 ^ 2^/2a, 

yif ~ 1 + aV2, yK+i 1 - aV2. (99) 



Note that if if = 1, then (3 > 2 and the above discussion is still valid, but we just have the two classes of solutions yi 
and y-jyi. 
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FIG. 2. Examples of a graphic resolution of equation (^) in the case f3 > 0. We have chosen — (0.01,0.1,0.5). The left 
plot describes the situation (3 — 2.3 and the right plot the case where /3 = 0.3. In the latter case K = 4. 



2. Oscillation probability 



Assuming that a/ (3 ^ 1, we can expand (73) to get the following forms for the coefficients /J.. Assuming that 
/3x - K 0, 



f 



f 

.1 X 



(1 - pj) 



c? 



\ (|/3k - 1| « 2\/2a) 

(^^^2 (|;3^ _ i| » 2V2a) 

i (|/3k-1|«2V2«) 
\-^T{i3-^) (|/3K-l|»2\/2a) 



(100) 
(101) 

(102) 
(103) 



When /3k - 1 > 0, is then given by ( |102D and /^^^^ is of the form (|103|) 
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It can be checked that, at lowest order ^f^. = 1- One can further check that when \(3k — 1| ^ 2\/2a the first 
correction to (101-102) is — a^I(/3~^)/2/3^ and that J2 fx =1 ^-l^o satisfied (which indeed has to be order by 
order). Now, the oscillation probability (p^) reduces to 



^^(7 ^ 5) 



A 1/2 



(1 - e)s\T? 



(|/3i^-l|«2V2a) 
{\Pk - 1| » 2V2a), 



(104) 



with rj = a^X(/3~^)//3^ as announced in (|63|). Thus, when \I3k — 1| ^ a, we are in a weak mixing regime and the 
probability is obtained by summing over all the individual probabilities. Otherwise, we are in a regime of strong 
mixing with the state and the oscillation length with this state is now given by 



A 1/2 ■ 



C. Summary and discussion 

In the limit where o;'^ <C 1, we have estimated the oscillation probability to be 



. /3 < 0: 



P(7->g)^4a2^- 



1- A 



• /3 > 0: 



^(7 ^ 9) 



/sin 

„2 



A 1/2 



i^K (l-/3.)'-i 



1-ft A 



4a' E^>l sin 



l-/3i 



A 



AW 



(|/3K-l|«2V2a) 
(|/?A— l|»2V2a2), 



with 1=1 — a^X(/3~^)//3^, as announced in (61) and (|6^). When |/3j — 1| ^ 2\f2a^ for all j, it is interesting to study 
the limit where u ^ £osc in which case we can assume that the sines can be replaced by their average value to get the 
two limiting behaviours: 



P(7 -> 5) ~ 2a2 1 



P(7 ^ 5) ~ TT 



45/32 



m » 1), 
«i). 



(105) 
(106) 



As explained in § [V, we respectively see on ( 105| -106) the small and large radius regimes where the extra-dimensions 
either have no effect ( |105| ) or enhance (106) the probability. We also find the regime of strong mixing with the K^^ 
KK graviton and the effect on the oscillation length as discussed below equation (|6^). 



VI. ESTIMATION OF THE PROBABILITY IN A SIX DIMENSIONAL SPACETIME 



Let us now turn to the physically more interesting case of n = 2 extra-dimensions. Now, equation ( p9|) cannot 
be solved exactly in general, but its solutions can be well approximated when the coupling between the photon and 
the graviton (this coupling is measured by Am in the matrix ( Pq) ) is small enough compared to the typical mass 
parameters of the mixing particles [i.e. the diagonal terms in (|26|)]. We follow the same lines as in the previous 
section. 
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A. 1 > ^ > 



1. Eigenvalues 



In that case, the plasma effects dominate over the vacuum polarisation in Ax; all the (3i are positive and one can 
easily show that the roots of (p9) are such that 



(107) 



We introduce k{i) the index such that P^^-^ is the closest Pi to yi. From ( 107 ) one has that k{i) d {i — l,i}- The 
eigenvalue equation for the root yt can be rewritten as 



yi-1 = a 



2 Sfc 



J^k{yi), 



where Tk is defined as 



Tu{y)^a^Y. 



(108) 



(109) 



T with no subscript denotes the function defined by the sum ( |109| ) taken over all indices j from j — Ito j — Nd . To 
finish, we introduce the index K such that Pk is the closest Pi to unity and then. 



(110) 



As in the five dimensional case, the determination of the roots yi has to be split in the three following cases: 
• i < K — 1: We first show that k{i) = i. For that purpose, we consider the function Ti.{y) defined by 

H{y) = l-y + T{y). (Ill) 

This function is strictly decreasing on ]Pi-i, Pi[ {Ti. vanishes only once in this interval in y = yi). Showing that 
H ^ ft+ft-i ^ > 0, is then enough to prove that k{i) = i. Since one has 



l-^^i±^>^ and Vfc 
2 - 2 



P^ + P^■ 



> 







using (C19), one obtains 



9 1 

< Q^sup i^Q', V/5yj < QQ'^ for y < 1 and p < 1. 



(112) 



(113) 



(the constants Q and Q' are defined in equation ( C19| ) of appendix Comparing (112) and ( |113D one sees 

> and then that k{i) — i. 



that for a smaller enough than P (namely 2QQ'^ < 1), one has Ti. j 
We will assume in the following the slightly stronger constraint 

^,2 



(114) 



Now, we set yi = Pi — ei, with ei > 0. Equation (|108| ) can be rewritten as an equation for e^, with Ti = Ti{yi)^ 

^2 „,2 „ 



P^-l-T, /3, - 1 - 

the positive solution of which is given to leading order by 



0, 



(115) 
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when a and (3 verify the constraint (114) 0. The eigenvalues are then given at leading order by 



1-A 



(116) 



(117) 



• i> K + 2: Using a similar line of reasoning as in the previo us c ase, one can show that k(i) — i — \ and then 
that Di is given at dominant order by (for a and /3 verifying (114)) 



2 -^t-i 
/3.-1- 



(118) 



• i — K,K + \: We first estimate the root y^- We assume that 1 e [/3_fc, /3i<-+i[ (similar conclusions can be 
obtained when 1 G /?;<-[), we have 



2 2' 



(119) 



As in the previous case, this is enough to show that yK & [Pk-i + 2 ^ S /^ifi 0- We set y^: = Pk — and 
Tk = ^kIiIk) with > solution of 



the positive root of which is 



1-(3k+ Tk 



ex 



C?Sk 



= 0, 



(120) 



-1 + Wi 



4a2 



SK 



{/3k - 1 - W 



(121) 



For a smaller enough than /3 one can consider the two limiting regimes (we will not consider here the 
intermediate case, in order to simplify the discussion) 



f a^sF if \Pk - 1| < 2a7sF 
" 1 if \f3K - 1| » 2ay^. 



(122) 



Let us now turn to the evaluation of the root yK+i- The discussion mimics the previous one. Assuming that 1 
is not too close to ^^^f^"^^ R, one can show that yx+i &]Pk,Pk + ^^+^~'^^" [ under the condition (114). Then 



considers (under the condition on a and /? of footnote W) the two limiting regimes 



we write yx+i ~ Pk — ^k+i, with ex+i < 0. eK+i is_solution of equation (120) with Tk = ^K{yK+i)- One 

(123) 



Pk + ol^/s^ for - 1| < 2q:^J^, 
^^+' ~ I 1 + -^(1) for - 1| > 2a^. 



to establish this we have use d (pld ) and (C21) and showed that (114) leads to |/3i — 1\ ^ Ti and |/3i — 1| ^ 2ay^ which 
in turn leads to the expression ( [ll6| ). In the rest of this section a 3> 6 means that a > 106 which we assumed to be enough to 
neglect b with resp ect a. 
^here we use (114) again. 

*One has here to impose a slightly stronger condition than the previous one: namely 4QQ'J^ < 1 in order to be able to 
neglect with respect to 2q^sx. This also insures that the expressions (122) are valid. 

/3jf+/3jc+l ' 



^namely 



J. 
10 • 



When this is not the case, one obtains the same results as in (125) when 



li > 2ay^. 
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2. Oscillation probability 
We now need to expand the coefRcients (M) in order to estimate the probabihty (|5 



• For i < K — 1, we have from equation ( p8\) and (117) 

(1-ft)' 



J Xi 



1 



with 



No 



Sk 



Using ( |C20| ), one can then show that under the condition ( |114| ) 

ii^^»max(l,ty,) 



(124) 



(125) 



(126) 



so that, at dominant order, 



(1-A)^ 



(127) 



• For i > _ftr + 2, we find in a similar way (and under the same condition) 

2 Si-l 



(i-A-i)^' 



(128) 



• For i — K,K + \, we distinguish the two regimes \I3k — 1| ^ 2a^/sK and \I3k — 1| ^ 2as/JK. 
Assuming that \(3k — 1| ^ 2q^sk, we find |^ the dominant contribution to f^^ and fxK+i to be 

2 SK 



f 

J Xk 



and ff. ~ 1. 



(129) 



When \I3k — 1| ^ 2ay^s/f , the dominant contribution to /J^ and fx^^^-^ are 



(130) 



Now, inserting these results in (|57|), we find that the oscillation probabihty is given to dominant order, for \Pk — 1| ^ 
2a^/sF, by 



2 sin2 



Xi - XK+l 



i=Nd 



4 y aV 
^ 1 



t (1-Ay 



■ sm 



1-ft 



A 



(131) 



This expression, as announced, is analogous to (^) and corresponds to the case when no KK state mixes strongly 
with the photon. 

Now, when \(3k — 1| <C 2a^/sK, one has 



^Using (|C2C|) , we show that this expansion is valid under the condition of the footnote 
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XK - Xk+1 



2 E E/. 



2 sin^ 



a^i - xp 



which can be rewritten as 



P(7 5) ~ (1 - 7^) sin2[AA/Vi:^w] + 4 ^ a' (1-/3 )^ 



(1 - A) 



A 



(132) 



(133) 



with 7] = J2i^K Q^^'5i/(1 — A)^ 1- As in the five dimensional case it was obtained by imposing that /^i = 1 order 
by order. This corresponds to the case where one KK state mixes strongly with the photon and again the oscillation 
probability is found to be equivalent to (|6l|). 



-1 < /3 < 



We now consider the case where f3 < (i.e. when the vacuum contribution dominates over the plasma effects in 
Aa). It is easy to see graphically (see figure |l|) that the Np) + 1 solutions of the eigenvalue equation ( p9| ) are such that 



2/1 > 1, Vi+i e ]/3j+i,/3j[, Vnd+i < Pnd- 



(134) 



1. Eigenvalues 



We do not detail the computation of the eigenvalues since it is similar to the former case. It is e ven simpler since 
now we do not have to single out the mode K (look for instance to the five dimensional case § V_A). 



We have only to assume the less stringent constraint than (114), lO^QQ' < 1, in order for the following expansions 
to be valid. Under this condition one can show that 



j^i^i + "'Er 



%>1 - 



1 - Pj^l 



(135) 



2. Oscillation Probability 

The coefRcients ( p8|) are then given at leading order by 

f2 ~ 1 

Jxi — -^1 



2 Sj_i 



(l-/3,_l)2 



The oscillation probability (^) can be expanded as 



F(7^g)^4/J^^4sin2 



where we have neglected higher order terms. This can be rewritten as 



P(7 ^ 5) - 



■ sm 



1-/3, 



A>u 



(136) 



(137) 



(138) 



which is again analogous to (pSj). 
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> 1 



We introduce the new parameters a = Aj\//A„i < and Pi = Am'V^™ — -P(ri)' 7 = ^a/A™ and z = a;/A,„. The 
eigenvalue equation (A9) can be rewritten as 



Nd 



Nd 



z — = a 



(139) 



One sees easily graphically that this equation admits one negative root zi and that the other Zi (i > 3) verify 
Zi &](3i-i, f3i[. We discuss here only the case where 7 is closer to than to 1 (we assume that 7 < 1/4). The discussion 
is very similar to the previous cases. Under the condition lOa^QQ' < 1, one finds that the roots Zi with i > 2 are 
given by 



a Si-i 



Zi^^ 13.1-1+ 5" "' \ so that fl^ 



Under the condition AaQQ' < 1, the two roots zi and Z2 are given, for 7 > by[^ by 



if 7 ^ 2\a\ and by 



from which we deduce that the coefficients are given either by 



if 7>2|a| 



4^f if 7«2|a| or by | {2' I 1^ if 7 » 2|a|. 



J X 



When 7 <^ 2|a|, the oscillation probability is given by 



X-2 — 7-^ 



P(7 ^ 5) ~ (1 ~ r)) sin2(AMM) + 4 ^ sin 



and when 7 ^ 2|a| by 



P(7-5)-4^- 



i>i 



■ sm 



A 



A" 



The small coefficient fj = J2i>2 / Pi is obtained as in ( |132D 



(140) 



(141) 



(142) 



(143) 



(144) 



D. Summary 

In the limit where < I we have estimated the oscillation probability ( p7| ) for a six dimensional spacetime to be 
• 1 > /3 > : For \I3k - 1| > 2a^JF, 

i=No 



P(7^g)c^4 ^ a'- 



t (1-A)^ 



■ sm 



1-A 



and for |/3a' — 1| <C 2ay/sK 



P(7 ^ 5) ~ (1 - 77) sin2[AMVs^u] + 4 ^ a^—ll— sin 



A 



AM 



(145) 



(146) 



these results being valid as long as AQQ'^ < 1 and the coefficient 77 being defined in equation (133) 



''For 7 < the results are similar; one has only to exchange the expressions of zi and Z2 in the case I7I ^ 2|a| 
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• -1 < /? < 0: 



• 2 

■ sm 



1-/3, 



(147) 



valid if 10^ QQ' < 1. 



P(7 ^ 5) - (1 - »?) sin2(AMii) + 4^] ^ sin' 



2 



■ sm 



1-/3^ 



A 



xu 



(148) 



respectivel y for 7 <C 2|a| and for 7 2|a|, the result being valid if lOa^QQ' < 1 and the small coefficient fj is 
defined in (p3|) . 



VII. MIXING IN AN INHOMOGENEOUS FIELD 

In all the previous sections, we have assumed that the magnetic field was homogeneous. This is however a very crude 
approximation for most of the realistic physical systems. In this section we first extand our analysis to inhomogeneous 
magnetic fields and give some implications of the inhomogeneity of the external field. 



A. Computation of the oscillation probability 

Following Q], we rewrite the equation of evolution ( p5| ) as a Schrodinger equation 

where we set V = (A, G^^-', . . . , G*^^-'). The two matrices TCq and Tii are respectively defined by 



Tio{u) = uj + 







.(1) 



and 



Hi(u) 



/ Am 

Am 



V Am 



Am\ 
I 



(149) 



(150) 



(151) 



We assume that Tii is a perturbation compared to Tio- This approximation is equivalent to saying that Aj\//AA and 
Aj\//A„i are small compared to unity, i.e. that a and a/ (3 <C 1. When Hq is inhomogeneous only Am and A>, 

depend on u while all the A^^ are constant. 



We first solve (149) at zeroth order, i.e. by neglecting Tii with respect to Tio, as 

V(")(w) ^ L/(w)V(°)(0), 
where the evolution operator U is defined by 



(152) 
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U{u) = ex.p—i / Ti.a{u')du' 



(153) 



Note that at this order there is no mixing effect since Hq is diagonal. 

Th e general solution of ( |149| ) is obtained by shifting to the "interaction representation" where Vint = C/^V so that 
(149) can be rewritten as 



■^intVi, 



int — ' Lint ''int 



with Tiint = U^TiiU. This equation can be solved iteratively by setting Vint ~ X^'^^iift with 

/ d^.'7^int(u')VrAu'), 



V 



(154) 



(155) 



with the initial condition v'>^2 = '^i^)- 

Introducing the basis {A^Qq} with A = (1,0, ...,0) and Qq being the state of the q^^ graviton, Q]^ — (5* for 
i e {1, . . . , + 2} and starting with an initial state describing a pure photon, i.e. V(0) = A{Q)A, we obtain 



-^^dz Am (z) V e^ "^^ ^(0)^; , 

Jo 



(156) 



where we have used that HqA = A^^, HiA — Am J2q Qq ^^"^ T^iQq ~ A^nGq- If we restrict to the first iteration, 
the oscillation probability is then given by 



^ 9) =Y.\{^^ v(°)+vw) 



V 



(1) 



(157) 



that is, using ( |l5q ). 



n Jo 



(158) 



We can check that in a homogeneous field we recover (63), i.e. the oscillation probability in the weak mixing case. 
Note that in the particular case of the weak mixing this method of computing the oscillation probability is shorter 
than the one used in the two former sections since it does not involve the determination of the eigenvalues of M. But, 
one has to assume that the probability is small compared to unity ji) , which is not necessarily the case for instance 
when we are in the strong mixing regime. 



B. Example of applications 



As an example, we consider the mixing in a periodic magnetic field of the form cos A^u with Aq > for which 
the oscillation probability, in the weak mixing regime, is given by (158 ) 



^(7->5)-E' 



dz Am cos 



3(Aoz)e " e -lo -"^ ' 



(159) 



Assume that Apiasma dominates so that we can neglect the variation of A\ with z then, the probability becomes 



1 



,(A^;;'^-Ai-^)^ 

where we have kept only the resonant term, which depends on the sign of Am''' — Ax, and where we have defined 

Ai±) = AA±Ao. (161) 



1 



{At'>-A[+^y 



^(n) _ ^(+) 



(160) 
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Now if lAm'-* — A^^^l <C |Am| or lAm'-* — ^| ^ \^m\ we find a strong mixing regime, meaning that because of the 
resonance there will exist a mode for which the probability is enhanced. 

In conclusion, the important scale that fixes the photon effective mass is now A^^-* ~ Ao if we are in a regime where 
|Aa| <C Aq. Wc can then have the same discussion as in the previous sections but with f3 defined as 

/3=^, (162) 

according to the sign of Aa. The length scale is now given by (loAo)^^^^ and it follows that we expect the two 
following effects: 

I. by increasing Aq we can hope to make /3 as small as wanted and thus to get a large enhancement of the 
oscillation probability. What happens is that the scale A^ is replaced by Ao and thus that a departure from 
the four dimensional case will be observ ed if (u>Ao)~^^^ < R. When the field is homogeneous, the scale A^^ is 



usually very large compared to R (see § VIII and § p^X| ) , which implies that there is little hope to see any effect 
of the extra-dimensions. By using an inhomogeneous field, we change the scale associated with the photon 
effective mass wich is now gouverned by Aq"^ that can be tried to be lowered to a scale close to R. 

2. whatever the sign of Aa, we expect to have strong mixing occuring for all values of Aq such that 

|A(;')-(Aa±Ao)|«Am. 
By varying slowly Aq or cj, we expect to see a series of strong and weak mixing regimes. 
The amplitude of these two effects will be discussed in the last section of this article. 



VIII. APPLICATION TO ASTROPHYSICS AND COSMOLOGY 



Magnetic fields are observed in most astrophysical systems but the origin of galactic and cosmological magnetic 
fields is still unknown . A possibility is that these fields have a primordial origin since such a magnetic field can 
be generated in a number of early universe mechanisms |4^ ] such as in collisions of bubbles produced in a first order 
phase transition |Q or during an inflationary phase 

The efficiency of the photon-graviton and of the photon-axion mixing depends both on the value of the magnetic 
field and on the spatial extension of this field, Ac say. We study the order of magnitude of these mixings on the cosmic 
microwave background, on pulsars and magnetars. 

The required quantities for our discussion are A^, A^, Apiasma and Aqed respectively given by equations ( [27| ) 
and (|3^). It is useful! to rewrite these quantities numerically as 



^.4x10-- 

1 cm ^ \ I G 

= 2 X 10-16 ( ^ 



(graviton 



A_ 







1 cm 1 

A, 



vlG, 
-2.5 X lO^s 



' /PQ ^ 

,1010 GeVy 



1 cm 1 



(axion), 

(2.5 X 1015)4/" VlTeV; VleV^' ' 



- -3.6 X I0-" 



Up 



I cm^i 

^QED 

1 cm^i 



: 1.33 X 10^27 I 



leV^ 
,f cV^ 



I cm^3 J 

2 



(163) 



where we have used the facts that 1 eV 



I G ~ 1.95 X 10 ^ eV in the natural Lorentz-Heaviside units 



where a = e'^/Air = 1/137 and the expression of the extra-dimensions radius 



R 



(2.5 X 1015)'/" 10-12 



D 



ITev 



eV" 



(164) 



We now restrict to the case n = 2. 
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A. Cosmic microwave background 



It has been shown that the isotropy of the cosmic microwave background (CMB) puts a hmit on the present value 
of a spatially homogeneous magnetic field to (Bq/IG) < 6.8 x 10~^{flQh^)^/^ A comparable bound has also 

been obtained for spatially inhomogeneous magnetic fields We study the magnitude of the photon-graviton 

conversion on the two following examples: 

• Large scales: we assume that we have a homogeneous magnetic field on the scale of the Hubble radius with 

77o =i 6 X 10"^ G. (165) 

The CMB photons are observed as a black body with a temperature of 2.7 K so that we approximatively have 



photons of energy piy52 



w ~ 10"^ - lO^^eV. (166) 

The caracteristic size of the system is the size of the Hubble radius 

Ac = 3000/i"^Mpc~ lO^^cm (167) 

where h is the reduced Hubble parameter. We also estimate the electronic density today to be about (see e.g. 
[1) 

ne~10"^cm-3. (168) 

• Degree scales: we assume a homogeneous magnetic field on the size of the Hubble radius at the last scattering 
surface. Since the magnetic field scales like (scale factor)^ and the energy of the photon as (scale factor)"^, we 
assume a magnetic field of 

Fo ~ 6 X 10-3 G (169) 

and consider photons of energy 

w~10"2-leV (170) 

at a redshift of z ~ 1000. The characteric size of the system is given by the Hubble radius at decoupling, i.e. 

Ac = 3 X 10^3/1-1 cm, (171) 

and the electronic density at the time of decoupling is of order (see e.g. fs^ ) 

Ue ~ IQ-^cm-^ (172) 

The main idea is that, since photons are converted into either gravitons or axions, some anisotropics must be induced 
on the scale of homogeneity of the magnetic field, mainly because of the angular dependence of the conversion rate. 
The effect between a direction parallel and direction perpendicular to the magnetic field must not exceed the observed 
CMB temperature anisotropy. The anisotropy of the CMB temperature between the directions perpendicular to the 
magnetic field (where the effect of mixing is maximum) and parallel to it (where there is no mixing effect) is then of 
order 



F(7 ^ g). (173) 



AT _ AT AT 
— - — ^ t' 

Observationally, we have the constraint that 

AT 

— < 10-^ (174) 

From figure ||, we deduce that in both cases, |Aqed| ^ lApiasmal so that A^ ~ Apiasma and thus (3 > 0. In the two 
considered regimes we have: 
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1. Large angular scale: 



a e 6.6 X [lO-l^lO~"] 



(175) 



Thus, we are always in a regime where a <C 1, /? > and \(3\ ^ 1 thus we expect at most effects of order 
which are completely unobservable. 



2. Small angular scale: 



a e 6.6 X [10^1°, 10~s] 

/3:^l.lxlOi^('^V 
^ V 1 TeV / 



(176) 



We are always in a regime where ^ 1, /3 > and 
observable effect. 



^ 1 and, as in the previous case, there will be no 



From this results, with see that /? is always too large to have any enhancement of the probability. Moreover in 
both cases the oscillation length, ^osc, with the lightest KK mode (as well as the oscillation length with any massive 
KK mode) is much smaller than Ac, and the mixing angle with the graviton zero mode is very small. The effects 
are the same as in a standard four dimensional spacetime and thus negligible Note that, in theory, we should 

have included the expansion of the universe but this will not change the result drastically. A detailed study of the 
photon-graviton mixing in an expanding four dimensional spacetime can be found in Q| and a discussion of the effects 
of the inhomogencity of the field in |q]. 
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FIG. 3. Aqed (dash line) and |Apiasmai (solid line) for the CMB on large (up) and small (down) angular scales. We see that 
we always have Aqed Apiasma. 



B. Pulsars 

As proposed by many authors (see e.g. |Q,|l^,Q), axions could be produced in the interior of neutron stars in 
nucleon-nucleon collisions. This would constitute the main cooling mechanism for these stars and thus puts limit on 
the axion production flux and mass. Such a production of KK gravitons in higher dimensional theories also exist and 
can be used to put bounds on the mass scale Md 

As originally proposed by Morris jl^ (see also [^), this axion (and now the KK gravitons) flux may be detectable by 
the secondary photons produced through the mixing with these particles in the neutron star magnetosphere magnetic 
field. These photons have a typical energy of 



27 



10"* eV, 



(177) 



i.e. of order of the average value of the neutron star interior temperature (about 50 keV). The primary photons can be 
well approximated by a black body spectrum with a temperature of Tiqs — 1 keV typical for the surface temperature 
of such stars. 

The idea is to detect a distortion of the star spectrum due both to the secondary photons and to the oscillation of 
primary photons. The typical value of the magnetic field in the neutron star magnetosphere is 



Ha ~ 10^2 G 

on a characteristic size of the system is of order of the neutron star size 

K ^ 10 km. 



Indeed, one cannot neglect the effect of the magnetospheric plasma and we estimate its density 55 56|| as 

-1 



7 X 10" 



IG 



(178) 
(179) 

(180) 



where P is the period of the pulsar and will be assumed to be about 1 second in the following. 
According to figure ^, we deduce that Aqed ^ | Api 

asma I SO tlia-t /? <C S-Ild ^\ — ^QED- TllGIl, it folloWS 



a ~ 3 X 10 
/5| ~ 3 X 10 



asma | 
14 



-8 / Md \ 



Tev) 



(181) 



and we are always in a regime where a <C 1, /3 < and |/3| < 1 and where the characteristic size of the system is far 
larger than the oscillation length. We expec t Mp to be of order 1 — 100 TeV, so that wc will get an amplification of 
order 1 — 10^ but still unobservable. From (181) we see that, contrary to the microwave background, the dominant 
length scale of the system is Aq^q so that 



^.2.4xl0-f^)-^ 
R VleV/ \1G) \lTeV J 



(182) 



which is smaller than unity for the typical value of magnetic field and wavelength considered here. By going to 
higher frequencies and higher magnetic fields we may get a larger amplification and thus a larger effect of the extra- 
dimensions, this is mainly the reason why we will turn to magnetars in the following paragraph. 

To finish, let us note however that in very strong magnetic fields one must take into account the photon splitting 
Ipl,^ which will compete with the photon-graviton mixing. We do not discuss this effect here. 
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log(O)) 

FIG. 4. Variation of Apiasma (solid line) and Aqed (dash line) respectively in function of the magnetic field Ho for a photon 
of tj = 10'' eV(top) and in terms of the frequency for a field of Ho = 10^^ G (bottom). 
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C. Magnetars and Gamma— Ray Bursts 



Magnetars are pulsars with superstrong magnetic field such as SGR 1806-20 [^sjjs^] where Hq ~ 8 x 10^^ G, i.e. 
two orders of magnitude higher than for ordinary radio pulsars. This object is associated with soft gamma ray bursts 
of energy of order 1 keV — lOOkeV. Other examples are GB790305 [^0[ and IE1841-045 and such observations 
are supported by models where the gamma ray bursts are triggered by cracking of the neutron star crust due to the 
magnetic stress []62| , |63| . 

So we consider a system such that 



1012_io15g, a, 



10 km, 



and 



10^ 



10^ eV 



Assuming that the electronic density is well approximated by ( 180| )^ we deduce that we are in the regime | Api 
Aqed as long as 



(183) 
(184) 

asma 



■ to 
.TeV 



IG 



> 4 X 10^ 



so that we can deduce that the QED contribution always dominates in such object and then that /? < 0. On figure ||, 
we depict the variation of a to show that we always have a <C 1. 



Now, effects of the extra-dimensions will appear when X^/R < 1 and this quantity varies typically from 10 



'18 



to 



10~^^ for n = 2 assuming Md ~ 1 TeV (see equation ( |182D ). Note also that u)^^ is of the order of R and that most 
of the effect comes from the fact that Aqed becomes large compared with R~^. Such objects may be interesting to 
detect the effects of the extra-dimensions but more data and a better understanding of soft gamma-ray bursts are 
needed before drawing any conclusions. 
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FIG. 5. Variation of a with respect to the plioton frequency for a magnetic field varying between Hq — 10^^ G and 
Ho = 10'^ G. 



IX. LABORATORY EXPERIMENTS 



Many experiments searching for light particles like axions were set up (see e.j 
review). We can classify the methods in the two following categories: 



12Hl4M,M and ma for a recent 



The pulsars cited above have a period ranging from 4 to 10 seconds. 
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• The direct methods in which a flux of axions coming from some astrophysical source (Sun, supernovae...) is 
tried to be converted into photons through an external magnetic field. These transitions have been used to put 
bounds on astrophysical axion fluxes and coupling constant |l^,|l^. One could think to use the same kind of 
experiments to put constraints in the case of a mixing with a large number of KK states. Let us first discuss the 
case of KK gravitons. The energy flux into KK gravitons from astrophysical object cannot exceed the bound 
on the energy flux in axions in the usual four dimensional case since otherwise the cooling rates of these objects 
will be too high . Let us then assume that the efficiency of detection is maximum and that the experiment is 
designed to collect all the emitted particles. Since each graviton is coupled with a much lower coupling constant 
than the four dimensional axion coupling constant accessible to these kinds of experiments, we do not expect 
that such direct detection methods will be able to see any KK graviton coming from astrophysical sources. In 
other words, since each KK graviton is coupled at tree level only to the photon (and not to other gravitons), 
there is no effect of the large number of KK states in these experiments. 

• The indirect methods where one tries to detect the mixing of the photon through its effect on a photon beam 
in a magnetic field, both on its amplitude and polarisation. Now, the photon being coupled at tree level to 
all KK states, one expects a departure from the usual case. The effect on the polarisation of the beam comes 
from the fact that for axions only the x component of electromagnetic wave couples to the axions. For the 
gravitons, both polarisations evolve according to the same equation but, due to the QED and Cotton-Mouton 
birefringence, A+ 7^ Ay which implies a phase shift between them. 

In the next paragraphs, we focuse on polarisation experiments to detect the phase shift. We first compute this 
phase shift in a Z? dimensional spactime and discuss two kinds of experiments respectively in a static and periodic 
magnetic field. We must emphasize here that the magnitude of the mixing with axions depends on the free parameter 
/pq (in contrast with the mixing with gravitons which magnitude is fixed by the value of the Planck mass), so that 
these experiments may be able to put constraints on bulk axion models. 



A. Phase shift in D dimensions 



Let us go back to the four dimensional case for axion. Then, from (^), we deduce that, starting from an initial 
state (^(0), G{0) = 0), the two polarisations + and x evolve respectively as, omitting a global phase cou, 



A+(w) = e-"^+"A+(0), Ax(w) = (e-^'^'x^cos^iJ + e-^'^^sin^i?) Ax(0). 
Now, restricting to the weak mixing regime where i? ^ 1, we can expand the A' defined in ( ^7|) as 



A'x ~ Ax +i?'(Ax - A„), 



a;, ^ A, 



^?'(A> 



A„ 



Expanding ( |185| ) to second order and taking into account (186) leads to 

A+{u) = e-*'^+"A+(0), Ax (u) = [1 - + i?2(e< - 1)] e-"^x"Ax (0) 
with C = (Ax — A„i)u. We deduce that the two modes evolve relative to each other as 



= [l-il?2^ + l?2(c< - 1)] 



)n^x(0) 



A+{u) L - ^ A+(0)- 

The relative phase and amplitude of Ax with respect to Aj^ then evolve as 



Ax 




1 - 2192 sin^ ^0 




Ax 


A+ 






A+ 



(0), 0(7.) ~ [(A+ - Ax)u - ^^iC ~ sinC)] m- 



(185) 

(186) 
(187) 

(188) 
(189) 



When we neglect the mixing effect (i.e. i? = 0), there is a phase shift due to the QED and Cotton-Mouton birefringence. 
The extra phase shift due to the fact that only one polarisation of the axion is affected by the mixing has been used 
to design experiments to put constraints on the axion parameters (see e.g. fl |l3| , p^ 64 1). 

In the case of a graviton, the two polarisations are mixed in the same way, so that the same computation leads to 



A> 



A+ 



1 - 2 



sin^ sm 



2fC+ 

2 



A> 



A, 



(0), 



.)- [{A+-A^)u + dl (C+-sinC+)-i?'x (Cx -sinCx)] 0(0) 



(190) 
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(we now have to keep the index X on d and on (). In four dimension, A^, = so that the phase shift depends only 
on the QED and Cotton-Mouton parameters. Its amphtude is proportional to Am so that it is roughtly 10 orders of 
magnitude lower than for axions. 

Let us now compute the phase shift in a -D dimensional spacetime. We assume that we are in the weak mixing 
regime and apply the method of § VII. As seen on (187), we must compute the solution of ( |l49|) up to second order. 
Following the same lines as for the computation leading to ( [156D , we can show that V^^^^ is explicitely given by 



(191) 



From which we deduce that, starting from a pure photon state, the polarisation A of the photon evolves as 



1-2^/ dyAM[y) dzAM(z)e-<y 



(192) 



In the case of a homogeneous field, we can extract from (192) the relative phase of the polarisation x with respect to 
the polarisati on + for the case of KK gravitons and bulk axions. In the latter case, only the polarisation x evolves 
according to (192) whereas the polarisation + evolves according to (187) so that 



(l){u) = (A+ - Ax)M-a2^ 



(Ax - A(;'))^i - sin(Ax - At;-)) J 0(0) (axion) 



(193) 



(A+- Ax)m-^s, 



[(Ax - A(^'^)w-sin(Ax - A^;;')) 



i>i 



(A+ - A^' V - sin(A+ - A^' V 



(l-/3(+))2 

with the notations used before. We split this result in three parts as 

4> = 0QED + 0CM + 4>KK 



(graviton) 



(194) 



(195) 



where the two first terms are the phase shifts due to vacuum polarisation and the Cotton-Mouton effect and are 
obtained by setting a = in (193-194). The third term is the specific phase shift associated with the mixing of the 
photon with either bulk axions or Kaluza-Klein gravitons. 



B. Polarisation experiments 

We now discuss two kinds of experiments designed to detect the mixing induced phase shift. As typical parameters 
we take w ~ 2 eV for the laser beam and a magnetic field which might be as strong as Hq — 10^ G. With these values, 
we have (for n = 2) 



J^c. 4x10-20, (graviton) 2 x 10"" ( -^f^ ) , (axion) 

1cm ^ 1cm ^ \10^"GeV ' 

AqeD ^q^in-17 ^plasma ^ ^■^■^^c^-n( \ ^"i _ t ^ i n-A, f 



/PQ 



3x10-1^ ^lPi2£iI^~-i.3xlO-iM--^ , -^l!^^ ~ -1.5 X 10-M -1^ . (196) 



4 



1 cm 1 ' 1 cm 1 V 1 cm ^ / ' 1 cm ^ \ 1 Tev 

Then, the QED and plasma effects are of the same order of magnitude but, since (/"qed S> 0cm, we neglect the 
Cotton-Mouton effect and define the phase shift ratio as 

nKK^^. (197) 

0QED 

One hopes to be able to measure a TZkk of order 0.1 B. 
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1. Multiple path experiments 



The idea is to make a laser beam reflect between two mirrors distant of /. Since the mirror are transparent to the 
axions and gravitons, the phase shift after N paths wiU be N(j){l) so that, in the case of axions, TZkk is given by 



7^ 



KK 



A, - A> 
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sin(A> 
(Ax 



At^^)l 



The decrease of the amphtude of the photon beams due to the creation of axions is 



61 ~ ia^NS^ — 
4- 1 



■ sm 



Ax - A. 



(n) 



(198) 



(199) 



Then, an enhancement of both the phase shift and the variation of the beam are expected due to the sums over all 
states. With the experimental values specified above, we have 

hi^io^(JiR.)\ (200) 

R VlTeVy ' ^ ' 

which is larger than unity. Then, we think that experiments in homogeneous magnetic fields will not probe the extra- 
dimensions since it will require to work with very high magnetic fields. Note that when we span the electromagnetic 
spectrum from the infrared to the X-ray, Xj/ R varies in the range 

2 



^~4x (10«-102) 
R 



Mp 
ITeV 



(201) 



2. Effect of a periodic field 



As seen in 



VII B 



one can hope to enhance the mixing effect by using a periodic magnetic field. For that purpose 
we need the pulsation Aq to dominate over A^ and /? to be smaU compared to unity. W ith t he previous numerical 
values, the first condition rewrites as Aq > 10~^^cm~^ and will be satisfied easily. Using (164), the second condition 
gives for a six dimensional spacetime 



A(7^ < 2.5 X 10^ 
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ITeV 



f— ) 

VleV/ 



(202) 



As stressed before, we have a departure from the four dimensional behaviour only if A-y < R. Now, since the magnetic 
field varies on a scale Aq'^, is gouverned by Aq instead of A\. Then an effect will appear only if we manage to 



create a field that can vary on scales of the order of the centimeter. 

In § VII B we also quoted the possibility of having a series of strong mixing regimes, which is specific of the existence 
of extra-dimensions. This requires that |Am''' — (A^ ± Ao)| ^ Am^/sI and can be performed either by varying Aq 
or Lu. Let us assume that Aq is fixed, since Aq 3> Am we have a strong mixing regime for the pulsations defined by 



,(p) 



2^2 An 



with a width of 



that is 
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leV 
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ITeV 
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M 
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A 



M 



1 cm 
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(204) 
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For instance if we assume that R is of order of the millimeter and that we consider a field varying on the order of the 
meter, we get for axions that 
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X. CONCLUSION 



We have used the property that, in an external magnetic field, the KK gravitons and axions couple at tree level to 
photons to show that there exists a mixing between the photon and these particles. We have computed this mixing 
and compared our result to the photon-graviton and photon-axion mixings in four dimensions. The main difference 
comes from the fact that the mixing matrix is now infinite and that a photon couples to a large number of massive 
particles. We then have discussed the physical implications of these phenomena in a general D dimensional universe. 
This leads us to conclude that for most astrophysical objects the effect of photon-KK gravitons mixing will be unob- 
servable. 



The main points of this study are: 

• We describe how to deal with the mixing between the photon and a large number of light particles. This extends 
the former results to the case of D dimensional universes where a photon can mix with all the Kaluza-Klein 
gravitons and possibly with bulk axions. In (|5^-|59|) , we gave the exact expression for the oscillation probability 
and we then discussed its amplitudes first qualitatively and then in a five and in a six dimensional spacetime. 

• When < R, i.e. if the caracteristic length scale associated with the photon wavelength and effective mass is 
smaller than the radius of the extra-dimensions, there is a departure from the four dimensional effect. Otherwise, 
the first KK mode is too heavy to be excited and everything, in general, reduces to the four dimensional situation. 

• Two limiting regimes have been found: 

— A large radius regime where the two following behaviours can appear 

* In the weak mixing regime, we have shown that the oscillation probability can be obtained by sum- 
ming over the individual oscillation probabilities and is then enhanced by a factor of order /3~"/^ in 
comparison with the standard four dimensional case. In that case the solutions can be found either 
by solving the eigenvalues equation or by considering the equation of evolution as a Schrodinger equa- 
tion and solving it iteratively in the interaction picture. The latter method generalises to the case 
of inhomogeneous magnetic fields but only applies if the oscillation probability is small compared to 
unity. 

* In the strong mixing regime the photon mixes preferentially with a given KK modes, which is possible 
if the plasma effects dominate over the vacuum polarisation. In that case a complete transition is 
possible. We note that this effect is more likely to happen in a D dimensional context than in a four 
dimensional spacetime and point out a specific effect of the extra-dimensions on the oscillation length. 

— A small radius regime where i? <C A-y so that the spacing between to KK modes is very large compared to 
the photon characteristic length. In that case, the probability is generally dominated by the contribution 
of the first state corresponding to the lightest particle and we are back to the four dimensional case. A 
consequence of this is that we can have an observable signature of the extra-dimensions only if A^ < R. In 
most of the systems this cannot be achieved, mainly because is very small; the only favorable situation 
happens in strong magnetic fields such as in pulsars and magnetars and when we deal with a magnetic field 
which varies on a small enough typical scale. 

• We have shown that, in the case of graviton, the effect of this mixing although enhanced is too small to be 
observed on the cosmic microwave background and on astrophysical objects such as pulsars. However, we point 
out that the effects can be larger with bulk axions. 



We discussed laboratory experiments designed for the search of axions in the light of this new framework and 



we computed the phase shift (193-194) between the two polarisations of a photon entering a magnetic field. As 
for the probability, we show that the phase shift is enhanced by the existence of the KK modes. In a periodic 
field, we show that the effect of the extra-dimensions can be important and that there exists a series of strong 
and weak mixing regimes that may be observed. More work to derive bounds on /pq is however needed before 
drawing any conclusion. 

On a more technical level, we have shown how to diagonalise a general matrix of mixing. This kind of matrices 
appears in different situations in extra-dimensions physics and these results can be used in many problems and 
in particular for neutrino oscillations. 
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• To finish, let us stress some important comments. First, we have assumed the vahdity of the Euler-Heisenberg 
Lagrangian to describe the vacuum polarisation. One has also to be aware that in strong magnetic fields such as 
in pulsar magnetosphere, photon splitting will be in competition with the photon graviton oscillation. We 

have not compared the strength of these two effects but we expect the latter to be dominant in high magnetic 
fields, such as in the magnetosphere of magnetars. We have also concentrated on gravitational waves even if 
there is also production of scalar waves. They are thought to be negligible, at least in the four dimensional 
case Concerning the axions, the effects may be more important than for gravitons, depending on the value 
of the coupling /pq, but we did not reconsider the bounds on the axion parameters. Further work is needed 
in that direction. We also stress that in general one needs a precise determination of their mass spectrum to 
compute the coupling of each mode to the photon. Most of these results were obtained for n = 1 or n = 2 
extra-dimensions for which the results respectively do not, or only weakly, depend on the cut-off of the theory. 
We also stress that depending on the exact physical situation the number of KK modes with which the photon 
can oscillate coherently can be drastically limited in comparison with the number of accessible KK modes. 
These decoherence effect implies that the UV cut-off of the theory A/max is expected to be much higher than 
the physical cut-off. A consequence of this is that the results obtained in the cases n ^ 1 and n — 2 can be 
extended to n > 2 without depending on Mmax- 



ACKNOWLEDGMENTS 

We wish to thank I'Ecole dc Physique Thcorique des Houches where this work was initiated, P. Binctruy, E. Dudas, 
J.F. Glicenstein, R. Lehoucq, M. Lemoine, J. Mourad, O. Pene, P. Peter and J. Rich for discussions. 



M.E. Gertsenshtein, Sov. Phys. JETP 14 (1962) 84. 

Ya. B. Zel'dovich and I.D. Novikov, The structure and evolution of the universe, Relativistic Astrophysic Vol. 2 (Chicago 
University Press) (1983). 

S.L. Adler, Ann. Phys. (N.Y.) 67 (1971) 599. 

G.G. Raffelt and L. Stodolsky, Phys. Rev. D37 (1988) 1237. 

L. Landau and E. Lifschitz, Theorie des champs (Eds MIR, Moscou), Quatrieme edition (1989). 

J.C.R. Magueijo, Phys. Rev. D49 (1994) 671. 

P. Chen, Phys. Rev. Lett. 74 (1995) 634. 

A.N. Cillis and D. Harari, Phys. Rev. D54(1996) 4757. 

J.E. Kim, Phys. Rep. 150 (1987) 1. 

G.G. Raffelt, Phys. Rep. 198 (1990) 1. 

M. Turner, Phys. Rep. 197 (1990) 67. 

P. Sikivie, Phys. Rev. Lett. 51 (1983) 1415. 

P. Sikivie, Phys. Rev. E52 (1984) 695. 

D. E. Morris, Phys. Rev. D34 (1986) 843. 

K. van Bibber et al, Phys. Rev. D39 (1989) 2089. 
C. Hagmann et al. Phys. Rev. Lett. 80 (1998) 2043. 
J.E. Kim, Phys. Rev. D58 (1998) 055006. 
P. Sikivie, Phys. Lett. B432 (1998) 139. 

G.G. Raffelt, Proc. of Neutrino '98, Japan, Y. Suzuki and Y. Totsuka eds. (1998). 
N. Arkani-Hamed, S. Dimopoulos and G.Dvah, Phys. Lett. B429 (1998) 263. 
N. Arkani-Hamed, S. Dimopoulos and G.Dvah, Phys. Rev. D59 (1999) 086004. 

E. Witten, Nucl. Phys. B471 (1996) 135. 
J.D. Lykken, Phys. Rev. D54 (1996) 3693. 

K.R. Dienes, E. Dudas and T. Gherghetta, Phys. Lett. B436 (1998) 55; ibid Nucl. Phys. B537 (1999) 47. 
L Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B436 (1998) 257. 

L Antoniadis, Phys. Lett. B246 (1990) 377; I. Antoniadis, C. Mufioz and M. Quiros, Nucl. Phys. B397 (1993) 515; I. 
Antoniadis, K. Benakh, Phys. Le tt. B326 (1994) 69; L Antoniadis, K. Benakh and M. Quiros, Phys. Lett B331 (1994) 
313; A. Pomarol and M. Quiros, | tiep-ph/980626:| . 



G.F. Giudice, R. Rattazzi and J.D. Wells, Nucl. Phys. B544 (1999) 3. 
T. Han, J.D. Lykken and R-J. Zhang, Phys. Rev. D59 (1999) 105006. 



34 



[29] P. Binetruy, C. D efFayet and D. Langlois, [iep-th/9905012|; P. Binetruy, C. Deffayet, U. EUwanger and D. Langlois, 
|iep-th/991021£ ] and references therein. 

V. Barger, T. Han, C. Kao and R-J. Zhang, Phys.Lett. B461 (1999) 34. 
S. CuUen and M. Perelstein, Phys. Rev. L ett. 83 (1999) 26 8. 

S. Chang, S. Tazawa and M. Yamaguchi, | tiep-ph/9908515|; tbtd., [iep-ph/9909240|. 



K.R. Dienes, E. Dudas an d T. Gherghetta, | |hep-ph/9912455 |. 
E. Dudas and J. Mourad, [iep-th/9911019|. 

J. Rich, Phys. Rev D48 (1993) 4318; C.Y. Cardall and D.J.H. Chung |iep-ph/990429lj ; B. Kayser, Phys. Rev. D24 (1981) 
110; W. Grimus and D. Stockinger, Phys. Rev. D54 (1996) 3414. 
W. Heisenberg and H. Euler, Z. Phys 88 (1936) 714. 

C. Itzykson and J-B. Zuber, Quantum field theory, MacGraw-Hill, New York (1980). 

L.D. Landau and E.M. Lifschitz, Electrodynamics of continuous media, MIR (1960). 

American Institute of Physics Handbook, 3rd. ed., D.E. Gray Eds. McGraw-Hil l, New York (1972). 

N. Arkani-Hamed, S. Dimopoulos, G. Dvah and J-M. Russcl, |hep-ph/98 11448 . 



K.R. Dienes, E. Dudas and T Gherghetta, |tiep-ph/9811428 [. 

L. Randall and R. Sundrurn, [ iep-ph/9905221 1; ibid, |hep-th/9906064 |. 

LS. Gradstheyn and I.M. Rizhik, Table of integrals series and product Ed. Academic, N.Y. (1990). 

P.P. Kronberg, Rep. Prog. Phys. 57 (1994) 325. 

K. Enqvist, Int. J. Mod. Phys. D7 (1998) 331. 

T.W.B. Kibble and A. Vilenkin, Phys. Rev. D56 (1995) 679. 

M.S. Turner and L.M. Widrow, Phys. Rev. D37 (1988) 2743. 

J.D. Barrow, P.G. Ferreira and J. Silk, Phys. Rev. Lett. 78 (1997) 3610. 

D.Puy and P. Peter, Int. J. Mod. Phys. D7 (1998) 489. 

K. Subramanian and J.D. Barrow, Phys. Rev. Lett. 81 (1998) 3575. 

G.F. Smoot et al, Astrophys. J. 396 (1992) LI. 

J.C. Mather et al, Astrophys. J. 354 (1990) L37. 

P.J.E. Peebles, Principles of Physical Cosmology Princeton University Press (1993). 
N. Iwamoto, Phys. Rev. Lett. 53 (1984) 1198. 

L. Shapiro and S.A. Teukolsky, Black Holes, White Dwarfs, and Neutron Stars, Wiley, New- York (1983). 
P. Goldreich and W.H. Julian, Astrophys. J. 157 (1969) 869. 

S.L. Adler, J.N. Bahcall, C.G. Callan and M.N. Rosenbluth, Phys. Rev. Lett. 25 (1970) 1061. 
C. Kouveliotou et al. Nature 393 (1998) 235. 

Kouveliotou et al, to appear in Astrophys. J. Lett., [ iistro-ph/980914C ] 
Paczynski, Acta Astron. 42(3) (1992) 145. 



Vasisht and E. Gotthelf, Astrophys. J. 486 (1997) L129. 

Thompson and R. Duncan, Month. Not. R. Astron. Soc. 275 (1995) 255. 

Thompson and R. Duncan, Astrophys. J. 473 (1996) 322. 

Maiani, R. Petronzio and E. Zavattini, Phys. Lett. B175 (1986) 359. 



35 



APPENDIX A: DIAGONALISATION OF M IN THE GENERAL CASE 



The goal of this appendix is to compute the eigenvalues and eigenvector of the matrix M , to diagonalise it and to 
explain how to compute the probability of conversion of a photon into a graviton. 
We consider the (TV + 2) x (TV + 2) matrix defined by 



/ Aa Aa/ 
' (r- 





M 



Am a1? 



■■• ••• 

A^;?) 



Aa/ 






V Aa/ 



(Al) 



•• aW/ 



We restrict to a finite matrix since there is a cut-off in the theory as discuseed in § PVB| and all the notations are 
defined in § ||. 

Let us stress that the diagonalisation of this matrice is a purely technical point that appears often in extra-dimension 
physics (see e.g. 41|). 



1. Characteristic polynomial 

These notations being fixed, we compute the characteristic polynomial V{x) of the matrix M defined by 

V{x)=dct{M-xlN+2), (A2) 



where In+2 is the {N + 2) x (iV + 2) identity matrix. Developping (A_2) with respect to its first column leads to 



N 



q=0 



N 
9=0 



where Dq is the determinant of the comatrix of the element (q + 2, 1) given by 



N 

D,(:r) = (-I)^Aa/ n i^i'^-^) 



From (A3) and (A4), we deduce that the characteristic polynomial is given by 



N 
9=0 



N 



Ax-x-A\,Y, 



(a(«) - X 



(A3) 



(A4) 



(A5) 



2. Eigenvalues 

The characteristic eigenvalue equation Vix) — has A'^ + 2 real solutions since M.^ being a symetric matrix, is 
diagonalisable. To find all his solutions, we rewrite (Rq) as 



V{x) =A{x)B{x), 



(A6) 



with 



i=l 

Bix) = {A,-x)l[[At^-x)-Al,Y^ n '^{^ 



No 



i=l 



Nd Nd 



i=i £=i,e^i 



(A7) 



We have two kind of eigenvalues: 
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• Am''': they are solutions of A{x) = and are of order Si — 1 but are not solutions of B{x) = so that they are 
solutions of 'P{x) =0 with order — 1 and thus are eigenvalues of Ai of the same order. 

This gives us '^fj'iisi — 1) = + 1 — Nd eigenvalues of A4. 

• Xi. they are solutions of B{x) = and since B{x) is a polynomial of order Nj:, + 1 and since Ai is diagonalisable, 
we must have Njj + 1 such eigenvalues. To find them, we rewrite B{x) as 



B{x) 



Nd 



X - A 



2 \ " f£ 

(Air.) 



(A8) 



Let us stress that Xi ^ Am'"* since otherwise the cancellation occuring in the first factor is offset by a divergence 
in the second factor. It follows that the Xi are solutions of 



Aa - a: = A 



Nd 
2 \ " 



(a(^'' - x) 



(A9) 



This solution can be found numerically but we can find the main properties of these eigenvalues graphically 
from which we deduce that (A9) has Nd + 1 distinct solutions that we order as 



<i<ND + l 



Xi < . .. < Xn+1- 



(AlO) 



(r) 

In conclusion, we have found the + 2 eigenvalues of M which split in n eigenvalues Am each with multiplicity 
Si — 1 and in Nu + 1 distinct eigenvalues Xi. 



3. Eigenvectors 



To determine the eigenvectors V solution of 



we set 



so that (All) reduces to the system 



MV = xV, 



V = {v,uo, . . .,un) 



Axv + Am 

AmV + Alfuq = XUq 



(All) 
(A12) 

(A13) 
(A14) 



li X = Am^\ the eigenvectors generate a subspace of dimension — 1 a basis of which is given explicitely by 

p 



1 



v/(p+l)(p + 2) 



^G,.+, + (p+l)G, 



^=0 



< p < - 2, 



(A15) 



where {A, (Gg)o<g<Ar} is the initial orthonormal basis where we have written A4 in (Al). One can check that 
this is an orthonotmal family, i.e. that 
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• li X = Xi, for each eigenvalue we have a subspace of dimension 1 generated by the unit vector 

" I „ . , -2 \ A^''"'^^ A^) ^^^^^ 

in the basis {A, {Gq)o<q<N}- It is easy to show that they satisfy 

We have given the exphcit form of the N + 2 eigenvectors of Al. It is worthwile noting that the eigenstates Vr^+p 
mix the different KK modes together while letting the photon unaffected whereas the eigenstates Vx^ mix the photon 
with the + 1 KK gravitons. 



APPENDIX B: PROBABILITY OF OSCILLATION IN THE GENERAL CASE 



To compute the oscillation probability between a photon and gravitons in a constant magnetic field, we follow the 
method by Raffelt and Stodolsky jij and solve the equation of evolution (25) as 



V(m) ^ e 



V(0), 



where V = {A, G''^\ . . . , G'-^^}. We decompose this vector on the eigenvectors basis as 

No Si-2 n+1 
i—1 p—0 2—1 



(Bl) 



(B2) 



where /ii,p(0) and /i(0) are N + 2 coefficients. Injecting this decomposition in (Bl), we obtain 



V{z) 



_i=l p=0 



ri+p ■ 



(B3) 



The probability of a photon to be converted in KK gravitons is obtained by considering the initial state V(0) 
{1, 0, . . . , 0} and by computing 



N 



9=0 



(B4) 



Since only the modes associated with the eigenvalues Xi mix with the photon, we deduce that 



^(7 ^ 5) = 1 - ^^(7 ^ 7) = 1 - 



2 ^ixiU 



(B5) 



where the coefficients are 



N 



1+E 



H 



Xi - A 



(9) 



-1/2 



(B6) 
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APPENDIX C: UPPER BOUND ON AND \gj\ 



The goal of th is ap pendi x isto give an upper bound on the absolute value of the two functions J^j{y) and Qj{y) 
(defined in and (|cT^) (|cT7| ) ) when n = 2. These majorations are used in section K^flto determine the solution 

of the eigenvalues equation ( |10^ ) as well as the oscillation probability in a small coupling limit. We also give an 
upper bound on Si. We first give a bound on 



2— ^l 

»2 



y-P^ 



and 



{y - A) 



2 ■ 



(CI) 
(C2) 



We recall that the (3i are defined by Pi = Pi^P = pfP, where pi is a pair (n^, rm) of integers, and pi is defined by the 
second equality. We assume all along this discussion that /3 is positive and we order the Pi as Pi < Pi+i- We stress 
that Pi=0. 



1. < 2/ < Al < |3^2 

For i > ii, each pf belongs to a unique interval 

(Ki + kf < p, < {p,, + fc + 1)2, (C3) 
where k is an integer. Mk, the number of such p^, is bounded by four times the surface defined by 



Pt^+k< ^Jxl+xl<p^,+k + l (C4) 

in the real plan (xi,X2), since we have at most four pairs for each unit square cell. Thus, one has 

Mk<^^r{p^,+k + l). (C5) 
One can then easily obtain a majoration in term of (p^^ + k): 

Mk<q{p^,+k), (C6) 



with q = 167r. On the other hand, for each pi satisfying (|C3|), \1/{y — Pi)\ is lower than |l/(j/ — {pi^ + k)'^P)\, from 
which we get the upper bound on J'ii,i.2iy) 

\^^..^M\<^1^ 2 — VP- (C7) 
^^ k=o [VPiPti + k)) -y 



kmax is defined such that p^^ + fcmax < Pi2 < Pn + ^max + 1- The r.h.s. of (C7) is nothing else but a Riemann sum 
associated with the function f{x) = x/{x'^ — y). Since f{x) is decreasing for all x"^ > y, we have 



P J^/^i{p^^^l) ^^-y 2/3 V P{p^^-lY-y 

Let us emphasize that ( |C8| ) assumes implicetly that y < /3(pii — 1)^. Otherwise, the contribution of the Pi such that 
Pp\ < Pi < P{Pn + 1)^ the sum ( ^l| ) can be bounded by 

^^P..- (C9) 



/^Pii - y 

Since /3(pii — 1)^ < y < Ppf^, we deduce that 
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Pn<[^ + 1], (CIO) 



so that we get the majoration 

With similar arguments, on can show that, when y < (3{pi-^ — 1)^, \Qii,i2{y)\ is bounded by 



\yn.Ay)\ - p (^2 _ yy < 2/3 - 1)2 „ y pi^p^^ + fc_)2 _ y ) ■ (^12) 



and that otherwise 



\gnMy)\ < (^f3^^1y)2 + 1) + /3fe, + fcL)^-J- ^^'^^ 

The bounds ( |Cll| ) and ( |C13|) are vahd also in the case where y < (3{pi-^ ~ 1)^- 

2. < ft, <(3,2<y 

In that case, following the same line of reasoning, we obtain respectively for Ti^^i^ and Qi^^i^ and any y satisfying 
the above condition 

•^ii,»2U/ < ^"TT + ^l^ ' (C'l^ 



a (7 



3. bound on jTj and C/j 



We now give a bound on the expression defined by 

No 

^ 2/ -A 

Si_ 

^=l,^^J (y - A) 



TM^a^ 5^ (C16) 

No 



For J defined by 

ViT^J, |y-ft| > (C18) 

One has thus Vi 7^ J, |y — /3i| > /3/2, and in particular /Spj^^ — y > /3/2 and ?/ — (3p'^j_i > (3/2. Using equations 
as well as the value of pmax givenFI in (^l|), one obtains easily the following bound on \!Fj\ and P|. 



®Note that the energy cut-off pmax has to be regarded as a maximum one fixed by the underlying quantum regularisation of 
the theory. However, decoherence effects can reduce drastically the number of KK modes that have to be taken into account. 
Thus, the bounds on Q and Q' in a more realistic case are likely to be lower than the one given here. 

'^"Similar bounds would also apply to J^{y) and Q{y) when y is so that Vi G [1, No], \y — fti > 13/2. 
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2 

I-FjI < sup (y^,Q') (C19) 
|ej|<Q^sup ^^,2) (C20) 



where Q is a constant of order 10^ and Q' a constant of order 10. 

4. bound on s; 

In the real plane the euclidian distance between two different pairs of integers is bounded by 1, so that the number 
of pairs of integers on a given closed curve is always lower than the length of this curve. It is then easy to obtain the 
following bound on Si which represents then number of different pairs of integers on a circle of radius pi. 

s^ < 2ttp, = 2^y^. (C21) 

To finish, let us note that the properties of the series Si have been studied by Gauss around 1800, see e.g. 
|ittp : //mathworld. wolf ram, com/rn.html for details and references. 
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